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temperature of superconductors is revisited. Anderson's pairing condition between time- 
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may find that the phonon-mediated interaction decreases exponentially by localization. 
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Impurity scattering can affect the transition temperature of superconductors by changing 
the electron density of states and the phonon-mediated interaction. It was shown 1 that the 
impurity effect on the electron density of states is overestimated in Abrikosov and Gor'kov's 
(AG) theory 2 of impure superconductors. In their comment 3 on ref. 1, Abrikosov and 
Gor'kov argued that their theory should be reinterpreted in terms of the Eliashberg theory 4 in 
order to satisfy Anderson's theorem. 5 Their argument follows from Tsuneto's 6 application of 
the AG theory to the electron-phonon model. His result apparently agreed with Anderson's 
theorem. The overestimation didn't occur because of the absence of the cutoff problem in 
the Eliashberg theory. 

However, the impurity effect on the phonon-mediated interaction is missed in both 
theories. 1 ' 7 Note that Anderson's theorem does not explain the exponential decrease 1 ' 8 of 
the phonon-mediated interaction by localization because it is valid only to the first order in 
the impurity concentration. 1 ' 7 The purpose of this study is to construct a strong coupling 
theory which can take into account the impurity effect on the phonon-mediated interaction 
precisely. The change of Coulomb interactions and phonon properties caused by impurities 
are neglected. 9 

First, Anderson's pairing condition is found from the physical constraint of the Anoma- 
lous Green's function. 10 Incorporating the pairing condition into the self-consistency equa- 
tion, one finds the exponential decrease of the phonon-mediated interaction caused by the 
localized wavefunction. 

We follow the real space formalism of the strong coupling theory by Eilenberger and 
Ambegaokar. 11 (See also refs. 12-16.) The Hamiltonian for the electron-phonon interaction 
takes the form 



where ^(r) and 0(r) are the electron and phonon field operators, g is the coupling con- 
stant. The equations of the motion for the thermodynamic Green's functions G(ou n , r, r') 
and F + (uj n , r, r') are given 




(1) 
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{iu n + -^J 2 + ^( r ) + v) G (u n , r, r') = 5(r - r') + g 2 J dr £(^„, r, r )G(u n , r c , r') 

+ ^ 2 y rfr o 0(w„, r, r D )F + (cj n , r G , r'), (2) 

(-^n + ^V 2 + V(r) + fi)F + (uj n , r, r') = # 2 | dr E + (u; n , r, r )F + (u n , r , r') 

+ # 2 y dr o 0+(w n , r, r )G(u;„, r , r'), (3) 

where 

S(w„,r,r') =T^ J D(w n ,^,r,r , )GK/,r,r'), (4) 

n' 

0(<j n ,r,r') = T^D(w„,w n /,r,r')F + (w n /,r,r'). (5) 

n' 

V( r ) — J2i V 5(r — Hi) is the scattering potential of the impurities and uo n = (2n + l)nT . D 
is the phonon Green's function. 

It is usually assumed that the electron-phonon interaction is local, 4 ' 11 ' 12 (i.e.), 

(u n , w„«, r, r') = 5{r- r')X(u n , u n >), (6) 
SK,r,r') = 5(r-r / )SK,r), (7) 

(f)(u n ,r,r') = 5(r - r')0(u; n ,r). (8) 

The generalization to the non-local interactions is straightforward and will be published 
elsewhere. 

The normal-state Green's function Gjv(cj„, r, r') satisfies the equation 

(iu n + t^V 2 + V(r) + n)G N (u n , r, r') = 5(r - r') + g 2 J dr E(cj n , r, r )G N (u n , r , r'). (9) 

From Eqs. (3) and (9) the Anomalous Green's function F + (uj n , r, r'), near the transition 
temperature, can be rewritten in the form 

F + (u n ,r,r r ) = g 2 J dr G N (-u n ,r ,r)(j) + (uJ n ,r )G N (uJ n ,r ,r f ). (10) 



Accordingly, we obtain the self-consistency equation for + 
+ (w n ,r) = T^A(w n ,uv)F + (^V,r,r) 

n' 

= 9 2 TJ2 Kun,u n >) j dr G N (-u n <,r ,r)G N (u n ,,r ,r)(j) + (uj n ,,r ). (11) 

n' 

In the presence of the impurities, the normal-state Green's function Gn may be given by 

G^„,r,rO = £^#^. (12) 

where Z(u n ) is the renormalization factor and ip m (j) is the scattered eigenstate. 

The pair potential A*(u n ,r) is defined by A*(o; n ,r) = 4> + (uo n) r)/Z(uo n ). Therefore we 
find the self-consistency equation for the pair potential to be 11 

A* ( u n ,r)Z(u n ) = 

g 2 T^2\(u n ,uj n ,) / dT G N (-u n >,r ,T)G N (u n >,T ,T)A*(u n >,T )Z(u n >). (13) 

n' J 

Let's compare this with the self-consistency equation in the Gor'kov's formalism which is 
given: 

A*(r) = VTJ2 I G(-^,r ,r)G(^,r ,r)A*(r )rfr . (14) 

Eq. (13) has the additional frequency dependence in the kernel. Notice that the spatial 
parts of the kernels in both equations are the same. 

However, it was pointed out 10 that Eq. (14) allows extra pairings which violate the 
physical constraint of the Anomalous Green's function. If we substitute Eq. (12) (with 
Z = 1) into Eq. (14), there are extra pairings between m | and m' j as well as Anderson's 
pairing between m f and fh {. m denotes the time reversed partner of the scattered state 
m. m' I does not include rh {. Because of the same spatial parts of the kernel, Eq. (13) also 
gives rise to extra pairings. In the Bogoliubov-de Gennes equations, the situations are more 
clear. 17 The corresponding unitary transformation leads to the vacuum state where pairing 
occurs between the states which are the linear combination of the scattered states. A similar 
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problem was found in Gor'kov and Galitski's (GG) 18 solution for the d-state BCS theory. GG 
allowed a superposition of several distinct types of the off-diagonal- long-range-order. 19,20,21 
Now we check the physical constraint of the Anomalous Green's function. If we average 
over the impurity positions, the Anomalous Green's function should be a function of r — r' 
due to the recovery of homogeneity. Using the scattered states VVn(jfc)( r )> 

^(r) = + J2 TT^T-JE e-^]e^> r + • • • , (15) 

q k k+q i 

it is easy to show that 

^ T ( r )^H( r, r P = ^^'^ + K 2 /'(r, r' ) + ■ ■ ■ 

^/(r-r'), (16) 

and 

U^)U-^') imP = ^ (r ' r ' ] + W( r -r>) + ... 

= <?(r -r'). (17) 

— * 

—imp means an average over impurity positions R{. /, /', g, and g' can be readily calculated. 
Note that even the zeroth order term of extra pairings between m j and m'(^ fh) [ is not 
a function of r — r'. To calculate the impurity average of the Anomalous Green's function, 
we need to consider the following average: 



F+(u n ,r,r') imp ~ i); n ^r)i)* mll (v')i) m ^r )i) mll (v )(P+(r ) imp 



+ correction terms. (18) 



From Eqs. (16)-(18), it is straight forward to show that 



F+K, r, r>r P ^ F+(uj n , r - r')^, (19) 



A*K,rf m VA*K), (20) 
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in the presence of the extra pairing terms. Only Anderson's pairing between m f and 
m I is compatible with the homogeneity condition of the Anomalous Green's function after 
averaging out the impurity positions. 

Consequently, the self-consistency equation needs a proper pairing constraint derived 
from the Anomalous Green's function. The resulting equation is 

A* ( LU n ,r)Z(uj n ) = 

g 2 T H^n, w n ') / dr {G N (-u n >, r G , r)G N (u n >, r G , r)} p , p A*(u n >, r )Z(u n >), (21) 

n' 

where p.p. denotes the proper pairing constraint. For ordinary impurities, Anderson's pairing 
condition between the time-reversed partners is obtained. The gap parameter A*(cu n ,m) is 
given by 22 

A*(u n ,m) = J ^ m (r)V4(r)A> n ,r)dr. (22) 
Finally, we find a gap equation 

A*Km)=^A(w m w nl )^V mm -p — A ^ w '' m ) (23) 



where 

V mm , =g 2 J |^ m (r)|V m ,(r)| 2 dr. (24) 

As a rule, the same result may be obtained by the k-space formalism. 4 
On the other hand, Tsuneto 6 obtained the gap equation 

= J^W J dq J d€ J dU e»-^V)^ ' (M) 

where n = 1 + ^pj) 00 = an d t is the collision time. Comparing Eqs. (23) and 
(25), we find that Tsuneto's result misses the most important factor V mm /, which gives the 
change of the phonon-mediated interaction due to the impurities. Anderson's theorem also 
assumes V mrn > = g 2 x l. 1 It is clear that this factor is exponentially small for the localized 
states, since the overlap of two localized states should in fact be exponentially small. 1 Let's 

6 



consider the states m and m' localized at r m and r m /, respectively. The localization lengths 
are taken to be £ and 0- The overlap of these wavefunctions are roughly proportional 
to \ip m (r m + Vm ~^ m ' )| 2 multipied by volume of the overlap region. It may be estimated that 

rm-r m /| 

V mm , ~ e — f d , (26) 

where d is the dimensionality of the system. Weak localization correction to the phonon- 
mediated interaction will be published elsewhere. 

In conclusion, we have constructed a strong coupling theory of impure superconduc- 
tors. Anderson's pairing condition, obtained from the physical constraint of the Anomalous 
Green's function, is incorporated into the self-consistency equation. One finds that the 
phonon-mediated interaction decreases exponentially by localization. 

I am grateful to Professor A. W. Overhauser for valuable discussions. This investigation 
was motivated by the response of Professor K. J. Chang at KAIST to my seminar on ref. 10 
on Nov. 19, 1994. This work was supported by the National Science Foundation, Materials 
Theory Program. 
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